Scale perturbation in valence bond ground states 
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A simple and efficient method for calculating the ground state for a class of antiferromagnet 
systems is presented. It combines the valence bond structure of the ground state for this class 
of systems and real space renormalization group. As an example, the entire Ifaldane phase for a 
spin-1 chain is described, from the AKLT model through to the Heisenberg model, ending at the 
critical WZNW SU2{2) model. The picture that emerges from this new method leads us to consider 
the relationship between the hidden topological order that characterizes the Haldane phase and 
potential applications in quantum communication technology. 

PACS numbers: 



One of the main problems in condensed matter physics 
is the formulation of accurate descriptions of the macro- 
scopic properties of a many-body systems. In order to 
achieve this aim, it is crucial to understand how do the 
low energy degrees of freedom emerge from the micro- 
scopic dynamics. Much of our knowledge about the 
emergence of these properties comes from the works by 
KadanoffQ and Wilsonl2l based on the Renormalization 
Group (RG) ideas (see also0). 

In the last decades, special attention has been given 
to the analysis of low dimensional antiferromagnetic 
systems where quantum fluctuations play an essential 
role (see 0] and references therein). In this context 
a new structure, known as valence bond, has emerged. 
This structure describes exactly the ground state of 
the Affleck-Kcnnedy-Lieb-Tasaki (AKLT) modelj^], con- 
structed from a Heisenberg-like Hamiltonian. Its inher- 
ent quantum nature, however, has brought new perspec- 
tives to the field of low dimensional antiferromagnetic 
systems extending beyond the AKLT model. 

The purpose of this work is to introduce a new method 
to characterize the ground states of quantum systems 
belonging to this family of states. This method takes the 
valence bond structure as a building block for ground 
states and combines it with real space RG in a simple 
way. We will use the quantum phase characterized by 
the AKLT model, the Haldane phase to illustrate the 
ideas coming from this method. 

To achieve these goals, we proceed in two steps: first, 
we will redefine the notion of the valence bond to take 
into account longer correlations in the system. Second, 
we will make a scale perturbation on the state to improve 
the correlations at every scale. The analysis of a spin- 
1 chain with this method leads to a simple picture of 
the Haldane phase. We will use this picture to describe 
the non-local properties of this phase0, @| in terms of 
entanglement Q. 

An essential requirement for any method that aims to 
describe the ground state of a quantum system is the 
efficiency in terms of computational resources and time 
needed to store the information about the state as well 



as to extract this information. 

Valence bond states (VBS) in one dimensional systems 
have a very simple and efficient description within the 
matrix product states (MPS) formalism[li|,|ll|. The goal 
of this formalism is to describe the wave function so that 
the entire information is contained in a product of local 
tensors: 



|*)^^AM|si)A[s2]k2)---AK] 



(1) 



where the index n points to a particular site in the ID lat- 
tice, {|sn)} is a basis in the local Hilbert space Hn — C'*, 
and A[sn\ = p{sn)T^'^"'^ are the tensors describing the 
amount of correlations in the state with p{sn) some prob- 
ability amplitude and {T^''")} a basis in Mdxd- Sev- 
eral properties have been characterized ^3 for this set of 
states; among them, we have that 

1. This set is a dense convex subset of translational 
invariant states. 

2. A transfer matrix that allows us to extract any ex- 
pectation values in an efficient way can be easily 
defined. 

3. If a state is the exact ground state of a system, 
then, it minimizes the energy locally i.e. the en- 
ergy density is equal to the lowest eigenvalue of the 
interaction. 

These three properties are the basis of this work. The 
first one allows us to use this class of states to formulate 
a variational ansatz; the second allows us to extract infor- 
mation from the system by a simple product of matrices; 
the third one shows that these states usually describe the 
infrared limit (defined later) for a given quantum phase. 

Another concept related with the third property is the 
notion of a parent Hamiltonian [slll^. This is the Hamil- 
tonian for which a variational wave function is the exact 
ground state. The parent Hamiltonian can yield some 
insight into the symmetries and properties of the phase 
under study; however, extracting the entire spectrum of 
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eigenvectors is usually a hard task and therefore to apply 
any perturbation method around the parent Hamiltonian 
is also difficult. With the method we are introducing, we 
will see that this task can be solved in an efficient way 
by defining a non-local MPS. 

The first goal we would like to achieve, is to describe 
the main properties of a quantum state in the neighbour- 
hood around VBS, i.e., when the parent Hamiltonian is 
modified with an irrelevant perturbation such that the 
system remains in the same phase. To accomplish this, 
we redefine the notion of a MPS to take into account 
longer correlations that might appear as the Hamilto- 
nian changes. The ansatz that we propose is to write the 
state as follows: 

1$) =^---A[so,Sl]|si)A[si,S2]|s2)--- (2) 

{s^} 

where A[s„_i,s„] — p(s„_i, s„)T'^^"\ with a non-local 
amplitude of probability p(s„_i, s„) and {T'^'*"^} a basis 
in Mdxd- 

This redefinition with a non-local MPS does not spoil 
the properties of the local MPS. In fact, this non-local 
redefinition reduces to the standard one when the proba- 
bility amplitudes are strictly local. Also, one can define a 
transfer matrix E[sn-i,Sn] — A* [s„_i, s„] (g) A[s„_i, s„], 
so that the information about the state can still be ex- 
tracted in an efficient way. Moreover, the non-local am- 
plitude p(s„_i, s„) can be used to apply variational meth- 
ods. Due to the fact that the correlations are attached 
to the physical indexes, it is straightforward to introduce 
symmetry arguments that simplify the problem. The 
drawback in this construction is that it is limited to sys- 
tems with short correlations. If we try to describe models 
with longer correlations by introducing more sites in the 
non-local MPS, the number of parameters grows in an 
exponential way and the method breaks down. 

Nevertheless, there is a way to avoid this problem: the 
combination of the VBS description of the wave function 
with real space RG methods. This is the second result of 
this work. 

The non-local MPS allows us to describe short range 
correlations. We can then improve the long range corre- 
lations in an efficient way, by integrating over the short 
distance degrees of freedom. This is nothing but a RG 
transformation. We therefore need to define the effective 
degrees of freedom and the dynamics at every scale in 
the transformation. 

The effective degrees of freedom are extracted from the 
non-local MPS by a coarse graining of a small cluster of 
contiguous lattice sites. For instance, the Hilbert space 
spanned by two lattice sites is given by 

T^'''^P{S2^,S2^+1)T^'"+''^\S2^, S2^+l) 

S2i,S2i+l (3) 



where {\si) = T,s2i,s2i+i C^if"'''*''^ l'S2i,S2i+i)} is a basis 
in the effective Hilbert space. 

The effective Hamiltonian is given by a real space RG 
method called contractor renormalization (CORE) cre- 
ated by Morningstar and WeinsteinTs^. By definition, 
this method keeps the eigenvalues of the low energy sec- 
tor and produces an optimal truncation operator from 
the original Hilbert space to the effective one. 

The low energy eigenvalues, {en}, can be obtained by 
exact diagonalization of the Hamiltonian for several clus- 
ters in the coarse grained lattice. In this work, we will 
see that it is enough to solve the interaction between two 
clusters or equivalently to diagonalize a four body Hamil- 
tonian. This step is the hardest computational problem 
of this method. 

We first calculate the eigenvectors, {|?^)}, for the four 
body problem and the effective Hilbert space for the con- 
tiguous clusters. CORE then shows that the optimal 
truncation operator is obtained by a Gramm-Schmidt or- 
thogonalization of the eigenvectors of the Hamiltonian 
projected on the effective Hilbert space, starting from 
the ground state upward. In this way, a basis, {|0ri)}, 
is built in which the first vector overlaps with the lowest 
energy eigenvector and those above, the second one with 
the second lowest and those above and so on, i.e. 

|0„) = 5] A™|m). (4) 

m>n 

Finally, the effective Hamiltonian is written as follows: 

-ffcff = ^en|0n)(0«|- (5) 
n 

At this new scale, we can then proceed to get the best 
non-local MPS and obtain the new effective system for 
the next scale. Usually, two situations can occur after 
several steps in the renormalization. The Hamiltonian 
either flows to a point where the interaction can be solved 
locally and the correlation length in the effective lattice 
model goes to zero, or the system is self-similar at every 
scale, the correlation length diverges and the mass gap 
goes to zero: at this point, the system is said to be at 
the critical point. 

This method combining VBS structure and real space 
RG transformation can be summarized as follows: 

1. Solve the short distance problem with the non-local 
MPS. 

2. Extract the effective Hilbert space from the MPS. 

3. Compute exactly the eigenvalues of the long dis- 
tance problem (the hardest computational step). 

4. Extract the truncation operator (CORE method). 

5. Obtain the Hamiltonian for the next scale and re- 
scale the unit of distance and energy. 
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6. Iterate the procedure. 

To illustrate the performance of this method, we will 
analyze the Haldane phase for a spin-1 chain. Several 
characteristics make this phase important on its own. 
First, the Heisenberg model is contained in this phase; 
this model is the subject of the Haldane conjectureQ 
which states that systems with half integer and integer 
spin have completely different behaviour. The former are 
gapless systems while the latter generate a mass gap that 
cannot be explained with perturbative methods, as this 
is a pure quantum effect. Second, this phase is a mag- 
netically disordered phase; it describes spin liquids with 
a hidden topological structure that cannot be repre- 
sented with the usual two point correlators0,0|. 

The Hamiltonian that we will study is written as fol- 
lows: 



H — ^ SiSi+i — {^SiSi- 



(6) 



where the operators 3°" , with a G {x, y, z}, correspond to 
the spin-1 representation of SU (2) and is a real param- 
eter. Depending on the actual value of /3, the behaviour 
of the system can change drastically. Three important 
points are described by this Hamiltonian: 

• [3 — \. At this value the system is integrable. It 
is described by the critical Wess-Zumino-Novikov- 
Witten (WZNW) model on the SU{2) group at the 
level fc = 2. This model happens to be equivalent 
to three massless Majorana fermionsfl^. 

• /3 = 0. This point describes the Heisenberg model, 
the subject of the Haldane conjecture. 

• /3 = — i. This value will be the starting point of 
our analysis. It is the AKLT0 model which has a 
VBS as the ground state. 

The AKLT model has all the ingredients that allow 
us to perform the scale perturbation as previously dis- 
cussed. Using a basis in the spin-1 representation such 
that 5"^|6) = ie"''=|c), with {a,6,c} = {x,y,z} and e"^'^ 
is the Levi-Civita tensor, the ground state of the AKLT 
model can be described locally as 



\AKLT), = 



V3 



(7) 



where a" are the usual Pauli matrices. This description 
corresponds to a MPS where p{si) = and T^**'^ — a^^ . 
Although this model is not the infrared fixed point of 
this phase, it is exponentially close to it'ls'l and contains 
the long range and topological properties of the Haldane 
phase. 

In the next step, we will modify the AKLT Hamilto- 
nian to obtain the Heisenberg model and introduce the 



non-local probability amplitudes p{si,Si+i) to describe 
the state at every scale. A basic analysis of the symme- 
tries that characterize this phase simplifies the problem. 
SU(2) and parity invariance reduce the number of pa- 
rameters in the variational amplitude. In addition to 
these symmetries, the normalization of the state is also 
required. 

The ground state energy per site and the correlation 
length will be completely characterized in our analysis for 
different values of the parameter /3 G [— i , 1] . Both quan- 
tities can be efficiently extracted from the MPS formal- 
ism. In the first plot, we show the ground state energy for 
the Heisenberg model (/3 = 0) and for the WZNW model 
{j3 = 1). The energy of the first model converges to 
Eq ~ —1.435522(5) and the correlation length converges 
to ^ ~ 3.8(0). These values can be contrasted with previ- 
ous numerical estimations where the energy is around 
3% higher. 

We can check how good this bound is by comparing 
the results obtained with another integrable model in 
the phase, the critical point. The second plot in FigQ] 
shows that the energy density of the WZNW converges 
to Eo ~ -4.175150(1). In Fig|21 the behaviour of the in- 
verse of the correlation length is plotted as we approach 
the critical point /3c — 1.03 ± 0.04, at which the correla- 
tion length diverges. 
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Figure 1: Plots of the ground state energy density as short 
distance degrees of freedom are integrated. By scale, we mean 
the number of steps in the renormalization, where we take 
groups of two contiguous sites every step. The first plot 
corresponds to the Heisenberg model where the energy per 
site converges to Eq ~ —1.435522(5). The second plot corre- 
sponds to the WZNW model where the ground state energy 
per site converges to Eq ~ —4.175150(1). 



Due to the theoretical analysis developed in we 
know that the correlation length behaves as ^ ~ |/3c — 
/3|^^, with = 1, as we change the coupling constant 
close to the critical value (3c and ^ ^ N'^, with i' — 1, 
as we modify the infrared cutoff N = 2'''^^'°. Fig|31 shows 
the plots for the behaviour of the correlation length. The 
slopes in the fitted curves show the scaling of ^ as the cou- 
pling constant approaches its critical value with a critical 
exponent 9 ~ 1.02 ± 0.02 and the scaling as the infrared 
cutoff is modified with h' ~ 1.00 ± 0.01. 

Finally, we would like to comment about the possi- 
ble consequences of this analysis on the hidden topolog- 
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Figure 2: Plot of the inverse of the correlation length as a 
function of the coupling constant /3. The points correspond to 
the numerical data, the solid line to the fitted curve and the 
dashed lines to a 95% confidence interval for the predicted 
responses given a linear fit. The inverse of the correlation 
length goes to zero at the critical value of the coupling con- 
stant ~ 1.03 ± 0.04. 
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Figure 3; Plots of the logarithm of the correlation length in 
base two. In both cases, the points correspond to the numer- 
ical data, the solid line to the fitted curve and the dashed 
line to a 95% confidence interval for the predicted responses 
given a linear fit. The first plot shows the logarithm of the 
correlation length as a function of log2(/3c — l3). The slope of 
the curve corresponds to a critical exponent ~ 1.02 ± 0.02. 
The second plot shows the logarithm of the correlation length 
as a function of the scale. The slope of the curve, in this case, 
corresponds to a critical exponent v ~ 1.00 ± 0.01. 



ical order of the Haldane phase'Tl and the potential ap- 
plication of this phase in quantum communication 
Verstraete and coauthors 9J showed that the behaviour 
of the non-local string order parameter of den Nijs and 
RommelsejS] can be understood in terms of the special 
entanglement properties of the Haldane phase. Using the 
VBS structure of the AKLT model, they showed that it is 
possible to create a maximally entangled pair with local 
operations and local measurements at any distance. The 
work presented here, shows that VBS structure does not 
change during the renormalization procedure and there- 
fore any ground state in the Haldane phase can be used 
to obtain a maximally entangled state, at any distance, 
using local measurement. 

To sum up, we have seen how to combine valence bond 
structure and real space RG. This method allows us to 
have complete control over the wave function as well 
as over the dynamics of the system. For that reason, 
we were able to extract thermodynamical quantities like 
ground state energy density and correlation length in the 
entire Haldane phase. Furthermore, we have shown that 



it is even possible to study critical points in an efficient 
way. 

Recently, three very interesting works [lij about simu- 
lation of quantum systems appeared, bringing new in- 
sight to the field. These works discuss yet other ap- 
proaches to this problem. 
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valuable remarks by Colbeck, Diir and Rodriquez. 
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